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, Abstract. We consider a multigravity approach to spacetime foam. As an 

<^> application we give indications on the computation of the cosmological constant, 

considered as an eigenvalue of a Sturm-Liouville problem. A variational approach 
with Gaussian trial wave functionals is used as a method to study such a problem. 
' We approximate the equation to one loop in a Schwarzschild background and a zeta 

■ function regularization is involved to handle with divergences. The regularization is 

closely related to the subtraction procedure appearing in the computation of Casimir 
energy in a curved background. A renormalization procedure is introduced to remove 
the infinities together with a renormalization group equation. 
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qh Despite of its recent interest, the term " Multigravity" is not completely new. Indeed, 

^ | in the early seventies, some pioneering works appeared under the name "strong gravity" 
or u f-g theory" [lj as a tentative to describe a sector of hadronic physics where a 
massive spin-2 field (the f-meson with Planck mass Mf ~ 1 GeV) plays a dominant 
role. Multigravity coincides with " strong gravity" or "/-</ theory" when the number 
of gravitational fields involved is exactly 2 ( u bigravity"). In this paper, we would 
like to use the Multigravity idea as a model of space-time foam[2] to compute the 
cosmological constant. Such a computation has been done looking at the foam as a 
large iV composition of Schwarzschild wormholes^S]. Nevertheless, the Multigravity 
framework seems to be more appropriate for such a computation. We recall that 
there exists a very large discrepancy between the recent estimates on the cosmological 
constant, which approximately are of the order of 10~ 47 GeV 4 , while a crude estimate of 
the Zero Point Energy (ZPE) of some field of mass m with a cutoff at the Planck scale 
gives E Z pe ~ !0 71 GeV i with a difference of about 118 orders^]. Let us see how to use 
multigravity, to approach such a problem. To this purpose we begin with the following 
action involving N massless gravitons without matter fields[5] 

N N 

S = J2S [Qi] = £ y^r / dV=ft [R (gd - A,] , (1) 

i=i i=i 1 J 

where Aj and Gi are the cosmological constant and the related Newton constant 
corresponding i th universe, respectively. Generally speaking, the total action should 



Space Time Foam 



2 



be of the form 



N 



Stat = ^2 S fed + ^ Sint (9i,92, ■ ■ ■ , gn) ■ 



(2) 



When A — > 0,the N world are non-interacting. This will be the context we are going to 
examine. In this way, the action So describes a gas of gravitons. Consider for simplicity 
the case N — 1 and the related Einstein field equations 



G^ u + Kg^v = 0, (3) 

where G^ v is the Einstein tensor. If we introduce a time-like unit vector such that 
u ■ u = — 1, then 



G^ u uV 



Ac (4) 

This is simply the Hamiltonian constraint written in terms of equation of motion, which 
is classical. However, the discrepancy between the observed cosmological constant and 
the computed one is in its quantum version, that could be estimate by the expectation 
value (A c ). Since 
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where R is the scalar curvature in three dimensions, we can write 
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where the last expression stands for 

[^}Jj,d 3 xH^[g tJ } 
V JV[g l3 ]^[g i3 }^[ giJ ] 

and where we have integrated over the hypersurface E, divided by its volume and 
functionally integrated over quantum fluctuation with the help of some trial wave 
functionals. Note that Eq.([7]) can be derived starting with the Wheeler-De Witt equation 
(WDW) [6 J which represents invariance under time reparametrization. Eq. (J7J) represents 
the Sturm-Liouville problem associated with the cosmological constant. The related 
boundary conditions are dictated by the choice of the trial wavefunctionals which, in our 
case are of the Gaussian type. Different types of wavefunctionals correspond to different 
boundary conditions. Extracting the TT tensor contribution from Eq. (|7j) approximated 
to second order in perturbation of the spatial part of the metric into a background term, 
(jij, and a perturbation, hij, we get 



A 



^ / d 3 x^G^ kl 



(2k) K 
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x , x k j ki + (2k) 



K X , X )iakl 



Here G^ kl represents the inverse DeWitt metric and all indices run from one to three. 
The propagator K 1 - (x, x) iakl can be represented as 



K 1 ( *\ 
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(9) 
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where h[^ >± ( af ) are the eigenfunctions of A2, whose explicit expression for the massive 
case will be shown in the next section, r denotes a complete set of indices and A (r) are 
a set of variational parameters to be determined by the minimization of Eq.flBJ). The 
expectation value of is easily obtained by inserting the form of the propagator into 
Eq.flU) and minimizing with respect to the variational function A(r). Thus the total 
one loop energy density for TT tensors becomes 



8nG W 
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(10) 



The above expression makes sense only for u>f (r) > 0, where u>i are the eigenvalues of A 2 . 
If we fix our attention on some particular background, for example the Schwarzschild 
background, the spin 2 operator A2, simply becomes 

The further step is the evaluation of Eq. ffTO]) . Its contribution to the Spin- two 
operator for the Schwarzschild metric will be 

(A 2 h TT Y i := -A S (A")J + « (l - ^) (h^J + 2 (Fth^t . (11) 

Ag is the scalar curved Laplacian, whose form is 

2MG\ d 2 f2r-3MG\ d I? 



A _f_2MG\ d 2 + /2r-3MG\ d 
\ r J dr 2 \ r 2 J dr 

and R°l is the mixed Ricci tensor whose components are: 

This implies that the scalar curvature is traceless. We are therefore led to study the 
following eigenvalue equation 

(A 2 h TT )i = u 2 h) (14) 

where uo 2 is the eigenvalue of the corresponding equation. In doing so, we follow 
Regge and Wheeler in analyzing the equation as modes of definite frequency, 
angular momentum and parity \7\. In particular, our choice for the three-dimensional 
gravitational perturbation is represented by its even-parity form 

(h even )i (r, tf, 0) = diag [H (r) ,K{r),L (r)] Y lm (tf, 0) . (15) 

Defining reduced fields and passing to the proper geodesic distance from the throat of 
the bridge, the system ffl4l) becomes 

-£s + l -^ + mj(r)]f l (x)=ul l f 1 (x) 

(16) 

■£ + l -^ + m 2 (r)]f 2 (x)=u J y 2 (x) 

where we have defined r = r (x) and 

m i ( r ) = U\ (r) = m\ (r, M) — m 2 , (r, M) 

■ (17) 

m 2 ( r ) = U 2 (r) = ml (r, M) + ml (r, M) 
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I (r, M) — > when r -> oo or r -> 2MG and m| (r, M) = 3MG/r 3 . Note that, while 



(r) is constant in sign, m\ (r) is not. Indeed, for the critical value f = 5MG/2, 
m\ (f) = and in the range (2MG,5MG/2) for some values of m 2 g , m\ (f) can be 
negative. It is interesting therefore concentrate in this range, where m\ (r, M) vanishes 
when compared with m\ (r, M). So, in a first approximation we can write 

m\ (r) ~ — m\ (r , M) 

(18) 

m\ (r) ~ +m% (r , M) 

where we have defined a parameter r > IMG and 7«q (r ,M) = 3MG/rQ. The main 
reason for introducing a new parameter resides in the fluctuation of the horizon that 
forbids any kind of approach. It is now possible to explicitly evaluate Eq. flTOl) in terms 
of the effective mass. By adopting the W.K.B. method used by 't Hooft in the brick 
wall problem|8], we arrive at the following relevant expression 
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mf (r) ) 2 

i — 1,2, where we have used the zeta function regularization method to compute the 
energy densities pi and where we have introduced the additional mass parameter /x in 
order to restore the correct dimension for the regularized quantities. Such an arbitrary 
mass scale emerges unavoidably in any regularization scheme. The energy density is 
renormalized via the absorption of the divergent part (in the limit e — » 0) into the 
re-definition of the bare classical constant A 

A _> Ao + A ^ = a + JL (ml (r) + m\ (r)) . (20) 

To remove the dependence on the arbitrary mass scale /x, it is appropriate to use the 
renormalization group equation. Therefore we impose that [9 J 

1 <9An T (LI,) (1 rprp 

& \ = P-rP T eJf (/*, r) , (21) 



where p^Jj (/x, r) is the renormalized energy density. Solving it we find that the 
renormalized constant Ao should be treated as a running one in the sense that it varies 
provided that the scale /x is changing 

A (/x, r) = A (/xo, r) + — (m{ (r) + m\ (r)) In (22) 

l07T /Xo 

The final form for the cosmological constant is|10j 

Ap (jVo) (r , M) / (r , M) 

SvrG _ 128tt2 M ^ 4 J [Z6) 
which has a minimum for 

ml (r ,M)y^ _ 1 

4/xi " ^ ( ] 
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We can now discuss the multigravity gas. For each gravitational field introduce the 
following variables (N, N{)^ and choose the gauge N- k ^ = 0, (k — 1 . . . N w ). Define the 



following domain D 
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and assume the following assumption: 3 a covering of S s.t. 



fc=l 



(27) 



for k ^ j. Then the Eq.((7]) turns into 
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(28) 



Each Sfc has topology S 12 x R}. Therefore, the whole physical space £ containing the 
energy density appears depicted as in the following picture A further simplification 




Figure 1. The space £ composed by the non overlapping spaces 



comes from the assumption that the different Newton's constants are all equal. This 
leads to a model which is composed by N w copies of the same world [TO] and on each 
copy the procedure contained between Eq. ffTOj) and Eq.fl23l) has to be repeated. Thus, 
the final evaluation of the "global" cosmological constant can be written 

■nax {1 £ AAg} + . . . + m ax y-j^ AA<t> 



A 



8ttG 



(29) 
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where A? is the eigenvalue obtained evaluating Eq.(|28l) on each The computation 
of the max is taken on each disjoint Note that in any case, the maximum of 
corresponds to the minimum of the energy density computed on the related hypersurface. 
It is interesting also to note that the whole procedure can be applied even in case of a 
massive graviton [TT] with a term of the form [121 02] 

m 2 r 

which is a particular sub-case of the Pauli-Fierz term[14j. 
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